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Abstract 

Spin-dependent pomeron effects are analyzed for elastic pp scattering and 
calculations for spin-dependent differential cross sections, analyzing power and 
double-spin correlation parameters are carried out for the energy range of the 
Relativistic Heavy Ion Collider (RHIC) at BNL. In this energy range, 50 < 
\/s < 500 GeV, the structure of pomeron-proton coupling can be measured 
at RHIC with colliding polarized proton beams. 
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I. INTRODUCTION 



The spin structure of the pomeron is still an unresolved question in the diffractive scatter- 
ing of polarized particles. There have been many observations of spin effects at high energies 
and at fixed momentum transfers |l]|§; several attempts to extract the spin-flip amplitude 
from the experimental data show that the ratio of spin-flip to spin-nonflip amplitudes can 
be non-negligible and may be independent of energy |3],^[]. In all of these cases, the pomeron 
exchange is expected to contribute the observed spin effects at some level. The large ra- 
pidity events that are observed at CERN [F3] and DESY ||, and all diffractive and elastic 
high energy reactions are predominated by pomeron exchange, thus making the pomeron a 
popular field of study. Extensive polarized physics programs are proposed at HERA, RHIC 
and LHC (see e.g. 0HH) i n order to shed light on these and other spin effects in hadron 
reactions. 

It is generally believed, based on calculations of the simplest QCD diagrams, that the spin 
effects diminish as inverse power of center-of-mass energy, and that the pomeron exchange 
does not lead to appreciable spin effects in the diffraction region at super-high energies. 
Complete calculations of the full set of helicity scattering amplitudes in diffraction region 
cannot be carried out presently since they require extensive treatment of confinement and 
contributions from many diagrams. Semi-phenomenological models, however, have been 
developed with parameters which are expected to be fixed with the aid of data from exper- 
iments [0,0]. 

Vacuum t-channel amplitude is usually associated with two-gluon exchange in QCD ||12|| . 
The properties of the spinless pomeron were analyzed on the basis of a QCD model, by 
taking into account the non-perturbative properties of the theory |13|J14|| . We refer to this 
as the standard-pomeron model in this paper. 

Some models predict non-zero spin effects as s — > oo and \t\/s — > limit. In these studies, 
the spin-flip amplitudes which lead to weakly altered spin effects with increasing energy are 
connected with the structure of hadrons and their interactions at large distances |IT||TT|j . In 
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reference [10|, the spin-dependence of the pomeron term is constructed to model rotation of 



matter inside the proton. This approach is based on Chou and Yang's concept of hadronic 



current density |lo[ . The model developed in reference []TJ] considers the contribution of a 
sea quark-antiquark pair to hadron interactions at large distances. 

This picture can be related with the spin effects determined by higher-order a s contri- 
butions in the framework of PQCD. Really, it has been shown in the framework of QCD 
analysis at fixed momentum transfer that these corrections to the simple two-gluon exchange 



may lead to the spin-flip amplitude growing as s in the limit of s — > oo. The similar 



effects can be determined by the quark loops contributions |17H 

The high energy two-particle amplitude determined by pomeron exchange can be written 
in the form: 

T{s,t)=isP{s,t)Vl hlP ®V^ F . (1) 

Here P is a function caused by the pomeron, with a weak energy dependence ~ (lns) n ; 
and V^ hP are the pomeron-hadron vertices. The perturbative calculation of the pomeron 
coupling structure is rather difficult and the non-perturbative contributions are important 
for momentum transfers of a few (GeV/c) 2 . 

The calculation of Eq. (0) in the non-perturbative two-gluon exchange model fl3| and in 



the BFKL model [18[ shows that the pomeron couplings have the following simple form: 

v£f = fa* r- (2) 

In this case, the pomeron contribution leads to a weak energy dependence of the differential 
cross section with parallel and antiparallel spins, and their difference drops as inverse power 
of s, leading us to conclude that the spin effects are suppressed as a power of s. 

This situation changes dramatically when large-distance loop contributions are consid- 
ered which lead to a more complicated spin structure of the pomeron coupling. As mentioned 
above, these effects can be determined by the hadron wave function for the pomeron-hadron 



couplings, or by the gluon-loop corrections for the quark-pomeron coupling [19(]. As a result, 



spin asymmetries appear that have weak energy dependence as s — > oo. Additional spin-flip 
contributions to the quark-pomeron vertex may also have their origins in instantons, e.g. 

IB- 



In the framework of the perturbative QCD, the analyzing power of hadron-hadron scat- 
tering was shown to be of the order: 

An oc ma s /ypt 



where m is about a hadron mass |22|. Hence, one would expect a large analyzing power for 
moderate pf where the spin-flip amplitudes are expected to be important for the diffractive 
processes. 

In this paper, we examine the spin-dependent contribution of pomeron to the differential 
cross sections with parallel and antiparallel spins, their possible magnitudes and energy 
dependence. We also estimate the possible experimental precisions for these observables in 
the RHIC energy domain. 

II. THE MODEL AMPLITUDES 

We use the standard helicity representation for the hadron-hadron scattering amplitudes: 
/i =< + + |M| + + > and / 3 =<H — \M\ H — >, helicity nonflip amplitudes; f 2 =< 
+ + \M\-- > and / 4 =<+-\M\-+ >, double-flip amplitudes; and f 5 =< ++\M\ + - >, 
single-flip amplitude. We assume, as usual, that at high energies the double-flip amplitudes 
are small with respect to the spin-nonflip one, f'2{s,t) ~ f<i(s, t) <C fi(s,t) and that spin- 
nonflip amplitudes are approximately equal, f + (s,t) = fi(s,t) ~ f 3 (s,t). Consequently, 
the observables are determined by two amplitudes: f + (s,t) and f-(s,t) = f$(s,t). These 
customary assumptions are also made for the models developed in P0|JTI| . 

We use the below normalization for the differential cross section: 
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and the analyzing power and the double spin correlation parameters are: 

-87T 

s(s — Am 2 ) 

0"o A NN = — — — 2|/_| 2 . 

s(s — Am z ) 

The measured spin- dependent differential cross sections can be written in the form: 

(3) 



iV( TT )/£ = 


o-l 


: TT ) 


= (0 O /4)[1 + A N {P 1 + P 2 ) 


i COS0 


+ A^PxPa cos 2 


N{ ll )/C = 


o-l 


: u ) 


= (0- o /4)[i-a 7V (Pi+p 2 ; 


1 COS(/> 


+ A NN P 1 P 2 cos 2 


N( tl )/C = 


01 


: n ) 


= (a /4)[l + A N (P 1 + P 2 ) 


I COS0 


- A NN PiP 2 COS 2 


iV( iT )/£ = 


0( 


: iT ) 


= (a /A)[l-A N (Pi + P 2 ] 


I cos</> 


- AnnPiPi cos 2 






0-Q 


= [iV(^) + iV(^) + N(^) 


+ iV( 


iT )]/£- 



£ is the luminosity and <Jq is the normalized differential cross section. Pi and P 2 refer 
to the degree of beam polarizations for the first and second beams, respectively, and <ft is 
the azimuthal scattering angle. The arrows indicate the transverse spin orientations of the 
interacting protons. 

If we adapt the following notation: 

of) = [AT(TT) + N{ n ))/C, a{ 1 ) = [iV( u ) + N(^)]/C, 

o-(TT) = [AT(TT) + N(U)]/C, <x( n ) = [N( n ) + N(^)]/C, 

then, the analyzing power, A N , and the double-spin correlation parameter, A NN , can be 
extracted from the experimental measurements; 

A _ ^( T )-^) _ Atr' -2Zm(/*/ + ) 

* 0-H + 0-W a |/ + | 2 + 2|/_p' lJ 
4 = ff ( TT )- ff ( n ) =^1= 2I/-I 2 ^ 

™ ff (TT)+ ff (U) a |/ + |2 + 2|/_|2" lJ 
Hereafter, Act s and Ao- d refer to the single- and double-spin cross section differences. We 
follow the model developed in |TT] closely and extend it further for the calculation of spin- 
dependent differential cross sections. Pomeron-proton coupling Vp pF is primarily connected 
with the proton structure at large distances and the pomeron-proton coupling looks like: 
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V£ P (p, t) = m Pll A(t) + 7 „B(t), (6) 

where m is the proton's mass, p is the hadron momentum and t is the four-momentum- 
transfer square. r y ll B{€) is a standard pomeron coupling that determines the spin-nonflip 
amplitude. The term mp^A(t) is due to meson-cloud effects. This coupling leads to spin-flip 
at the pomeron vertex which does not vanish in the s — > oo limit. Using Eq.(^), we can 
calculate the spin-nonflip and spin-flip amplitudes from the pomeron-proton vertex: 

\f + (s,t)\<xs\B{t)\, 
\f„(s,t)\cxmJ\t\s\A{t)\. (7) 



Hence, Vp pF determine A N and A NN parameters. Both of the above amplitudes have the 
same energy dependence and the ratio of spin-flip to spin-nonflip amplitudes in this picture 
gives: 

m|/_(M)l m 2 \A(t)\ A AK , n n „ ^ Ul . ftE ,™^ 2 



_ _ , r , _ 0.05 to 0.07 for \t\ ~ 0.5 (GeV/c) 2 (8) 

which is consistent with other estimates |§. 

The amplitudes A and B have a phase shift caused by the soft pomeron rescattering. As 
a result, the analyzing power determined by the pomeron exchange 

A N ^f = 2m^\ Im{AB*) (9) 

and appears to have a weak energy dependence. 

The model also takes into account the contributions of the Regge terms to both /+ and 
/_ amplitudes. So, the scattering amplitudes are 

f ± (8,t) = i8\J£ + fa + C J^±A f r_e 9] = is[fP± + f±l (1Q) 

where f±, f± 9 are functions which weakly depend on energy and c; are parameters. The 
asymptotic terms f± and f± 9 were calculated in the framework of the model [11]. The 



Born amplitudes in the form of Eq.(^) are modified by pomeron rescattering. The Regge 



contributions, Eg. (|T0D , were represented in the simplest exponential form. These, and some 
of the asymptotic function parameters, were obtained from the fit to the experimental data 



(for details see [ITJ). The model quantitatively describes all the known experimental data of 
the proton-proton and proton-antiproton scattering, from y/s = 10 GeV up to \fs = 1.8 TeV 
TT| , |2"4|| . We thus expect the predictions for differential cross sections at RHIC energies to be 



reliable. We neglect a possible odderon contribution in these calculations. 



III. SPIN CORRELATION EFFECTS 



The meson-cloud |TT] and the rotating matter current [10] models quantitatively describe 
the experimental data on elastic pp scattering at fixed momentum transfers and can predict 
physical observables (cross sections and asymmetries) at high energies. The predictions, 
however, differ in size and sign for asymmetries. 

We calculate the spin-dependent cross sections using the above described amplitudes in 
Eq. (|i0|) , and we use the parameters of the RHIC beams || for the estimation of statistical 
errors. We assume at \/s ~ 200 GeV, the luminosity is £ = 3 x 10 31 cm -2 sec -1 and 
the average beam polarization per beam is 70%. The typical geometrical acceptance of 
the detector is taken to be 20% and the running time is about a month. The momentum 
transfer binning, we take At = 0.05 (GeV/c) 2 at \t\ < 1.3 (GeV/c) 2 and At = 0.1 (GeV/c) 2 
at \t\ > 1.3 (GeV/c) 2 . 

The energy dependence of the real-part of the nonflip amplitudes around diffraction 
minimum (Imf + (s,t) ~ 0), is model dependent and may lead to different polarization 
predictions in this momentum transfer range. 



The local dispersion relations have been used in |TT] to determine the real-part of T(s, t). 
The model amplitude obeys the s — u crossing symmetry which permits us to describe 
quantitativaly all the specific effects in the elastic proton-proton and proton-antiproton 
diffraction scattering in the wide energy region (9.8 GeV < y/s < 1800 GeV). For 
example, model discribe the difference of the differential cross sections for proton-proton and 
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proton-antiproton scattering in domain of the diffraction minimum, the known polarization 
phenomena at y/s = 9.8 GeV and y/s = 52.8 GeV and the experimental data of proton- 
antiproton scattering at y/s = 540 and 630 GeV. All there quantities are sensitive to the real 
part of the scattering amplitude. So, we can believe that Re[T(s,t)} determined correctly 
and practically model-independent. 

The differential cross section calculations are shown in Fig. 1 a,b. The estimated er- 
rors are statistical and are less than 1%. At y/s = 50 GeV, the model described here 
quantitatively reproduces the ISR data ||23|| . The diffraction minimum defined by the zero 
of the imaginary part of the spin-nonflip amplitude is filled by the contributions of the 
real-part of spin-nonflip and spin-flip amplitudes. At y/s = 120 GeV, the dip becomes a 
shoulder and increases by an order of magnitude when the energy goes from y/s = 120 GeV 
to 500 GeV. The model gives the same asymptotical predictions for the proton-proton 
and proton-antiproton differential cross sections. Hence the cross section prediction at 
y/s = 500 GeV approximately corresponds to the SppS data at y/s = 540 GeV. 

Fig. 2 a,b show calculations for Aa d as determined in Eq.(^). It can be seen that the 
shape of Aa d is similar to the spin-averaged differential cross sections. It has a sharp dip 
and the magnitude of Aa d at the dip grows by an order of magnitude in the energy interval 
120 < y/s < 500 GeV. Although the dip positions for differential cross sections and that of 
Aa d seem to coincide at y/s = 50 GeV but at higher energies they move apart from each 
other since the dip of Aa d moves more slowly. The position of the Aa d minimum strongly 
depends on the model parameters. However, in the ranges of \t\ ~ 0.7 to 1.0 (GeV/c) 2 and 
\t\ > 1.8 (GeV/c) 2 , the results weakly depend on these parameters. In the regions far from 
the dips, both cross sections change slowly, especially at \t\ > 2 (GeV/c) 2 . 

The energy dependence of f + (s,t) and amplitudes are determined in Eq.(^) for 

the V pp p vertex. Hence, we have the same energy dependence for Aa d and cr(^) for the 
spin-pomeron models, and the ratio of these quantities will be almost energy independent, 
except around the diffraction minimum. 

The calculated ratio of the spin-parallel and antiparallel cross sections, R(s, t) = 
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c( )j shows only a logarithmic dependence on energy in its first maximum and is 

almost energy independent at \t\ = 0.7 (GeV/c) 2 , see (Fig. 3 a,b). It is clear from these 
results that the spin effects can be sufficiently large at all RHIC energies. 

Relative errors in R(s, t) are shown in Table I at two selected momentum transfers and 
at t m where the first maximum of R(t) is observed. It is worthwhile to note that the errors 
with the growth of energy decrease due to increase in the differential cross sections. Also 
note that in all standard pomeron models this ratio is predicted to be 

R( s ,t) =cr( n )/(7( n ) -> 1 as s^oo. 

However, the experiment data show a large deviation from unity [|TJ. It is also interesting 
that at low-energies (j>l = 11.7 and 18 GeV/c) p5| , p6| , R(s,t) does not change much with 
energy at \t\ ~ 2 (GeV/c) 2 and is similar to our predictions of R(s, t) at the first maximum. 

Our calculation for An is shown in Fig. 4 a,b. The magnitude and the energy dependence 
of this parameter depend on the energy behavior of the zeros of the imaginary-part of spin- 
flip amplitude and the real-part of spin-nonflip amplitude. 

The maximum negative values of coincide closely with the diffraction minimum (see 
Figs. 1 and 4). We find that the contribution of the spin-flip to the differential cross sections 
is much less than the contribution of the spin-nonflip amplitude in the examined region of 
momentum transfers from these figures. A^ is determined in the domain of the diffraction 
dip by the ratio 

A N ~Imf_/Ref + . (11) 

The size of the analyzing power changes from —45% to —50% at y/s = 50 GeV up to 
—25% at yfs = 500 GeV. These numbers give the magnitude of the ratio Eq.(|TT|) that 
does not strongly depend on the phase between the spin-flip and spin-nonflip amplitudes. 
This picture implies that the diffraction minimum is filled mostly by the real-part of the 
spin-nonflip amplitude and that the imaginary-part of the spin-flip amplitude increases in 
this domain as well. 
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We observe that the dips are different in speed of displacements with energy from Figs. 
1 and 2. In Fig. 4, one sees that at larger momentum transfers, \t\ ~ 2 to 3 (GeV/c) 2 , the 
analyzing power depends on energy very weakly. 

The expected errors for the analyzing power are small in nearly all momentum transfer 
ranges examined. They are summarized in Table II. 

As pointed out before, the model 0] predicts a similar absolute value but opposite sign 



for An near the diffraction minimum. The future PP2PP experiment at RHIC |§ should be 
able to provide data and help resolve these issues on the mechanism of spin-effect generation 
at the pomeron-proton vertex. 

IV. CONCLUSION 

In the framework of the standard-pomeron model, the spin-flip amplitude is defined only 
by the secondary Regge poles and the ratio 

K TT ) - *( U )]M n ) <x 1/a 
that rapidly decreases with growing s due to the standard energy dependence of the spin- 



flip amplitude. If we drop the asymptotic term in the spin-flip amplitude from Eq. (TTTJ) and 
keep only the second Regge terms that fall as 1/y/s, we obtain the pre-asymptotic Regge 
contributions to the analyzing power and the difference of the polarized cross sections (the 
results are shown in Fig. 5 and 6). 

The spin-flip amplitude of the second Regge contributions has a large relative phase 
compared to the spin-nonflip amplitude of the pomeron. Under this condition, too, the 
analyzing power can be large. As one observes from Fig. 5, the spin-flip amplitude defined 
by the Regge contributions can describe the experimental data at yfs = 23.4 GeV and 
give large effects at y/s = 50 GeV. At higher energies, however, the effect quickly falls 
and becomes insignificant. Note that at y/s = 50 GeV, can be positive and that 
its magnitude greatly depends on the size of the asymptotic term which gives a negative 
contribution to the analyzing power. 
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If the spin-flip amplitude is determined by the standard-pomeron model, its contribution 
will be clearly seen in experiments (see Figs 5, 6) if performed. In this case, the minimum 
in Aa d is at small momentum transfers \t\ < 0.7 (GeV/c) 2 , and is quickly shifted with 
growing energy. In the ordinary dip region, there is a maximum in the difference of the 
cross sections which falls as the inverse power of energy. Moreover, at \t\ > 2 (GeV/c) 2 , 
in the spin-pomeron model, cr and Aa d do not appreciably change with energy at fixed t 
(Figs. 1,2). The difference spin- dependent cross sections fall with energy as quickly as the 
increasing momentum transfers (Fig. 6). 

The energy dependence of cross sections cr(^) and cr(^) can be studied experimentally 
at RHIC. Note that significant spin effects can have small relative errors at momentum 
transfers \t\ ~ 2 to 3(GeV/c) 2 and direct information about the nature of the spin-flip 
effects in the pomeron-proton coupling can be obtained. The future PP2PP experiment at 
RHIC should be able to measure the spin-dependent cross section with parallel <r( TT ) and 
antiparallel cr(^) beam polarizations in proton-proton scattering and the energy dependence 
of the spin-flip and spin-nonflip amplitudes can be studied in the energy range of RHIC, 
50 < v/i < 500 GeV. 

The spin-structure of the pomeron couplings are determined by the large-distance gluon- 
loop correction or by the effects of hadron wave function. Tests of the spin structure of QCD 
at large distances can be carried out in diffractive reactions in future polarized experiments 
at HERA, RHIC and LHC. 
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TABLES 

TABLE I. The expected relative errors for R(s,t) = a(^)/a(^) are tabulated below at 
ti| = 0.5 (GeV/c) 2 , |*2| = 2.5 (GeV/c) 2 and at t m , where R(s,t) is maximum. 



v^[GeV] SR(h) SR(t m ) SR(t 2 ) 

50 0.05% 5.5% 4.1% 

120 0.07% 2.2% 4.1% 

250 0.07% 1.1% 3.7% 

500 0.09% 0.6% 3.6% 



TABLE II. The expected relative errors for An are calculated at |£i| = 0.5 (GeV/c) 2 , 

I ^2 1 = 2.5 (GeV/c) 2 and at t m , where An is maximum. 

v^[GeV] 5A N (h) 5A N (t m ) SA N {t 2 ) 

50 3.1% 3.7% 16% 

120 2.3% 2.0% 20% 

250 3.5% 1.4% 18% 

500 12.4% 1.1% 17% 
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FIGURES 

FIG. 1. a) The calculated differential cross sections of the pp elastic scattering at 
s 1 / 2 = 50 GeV (solid curve) and at s 1 / 2 = 120 GeV (dashed curve); b) s 1 / 2 = 250 GeV (dot- 
ted curve) and at s 1 / 2 = 500 GeV (dashed curve). The energy range between 50 GeV and 500 
GeV correspond to the entire RHIC energy for polarized protons. 

FIG. 2. a, b) The calculated difference of the double spin differential cross sections at the 
RHIC energy range. The error bars indicate possible statistical errors for a realistic experiment at 
RHIC. 

FIG. 3. a, b) The calculated ratio of the cross sections with parallel and antiparallel spins at 
the RHIC energy range. The error bars indicate possible statistical errors for a realistic experiment 
at RHIC. 

FIG. 4. The calculated analyzing power of the pp elastic cross sections for the weak energy de- 
pendence obtained in the model for the ratio |F +_ (s, 0)|/|F ++ (s, 0)| oc /(Ins) a) at s 1 / 2 = 50 GeV 
(solid curve) and at s 1 / 2 = 120 GeV (dashed curve); b) s 1 / 2 = 250 GeV (dotted curve) and at 
s 1 / 2 = 500 GeV (dashed curve). The error bars indicate expected statistical errors for a realistic 
experiment at RHIC. 

FIG. 5. The analyzing power for the pp elastic cross sections for the rapid energy depen- 
dence of the ratio |F +_ (s, 0)\/\F ++ (s, 0)| oc s' 1 / 2 (the experimental data at P L = 300 GeV/c 
are represented by open circles) ; the results of calculations are shown for Pl = 300 GeV (solid 
curve); at s 1 / 2 = 50 GeV (long-dashed curve); at s 1 / 2 = 120 GeV (short-dashed curve); and at 
s 1 / 2 = 500 GeV (dotted curve). 

FIG. 6. The calculated difference of the double spin differential for the rapid energy depen- 
dence of the ratio |F + ~(s, 0)|/|F ++ (s, 0)| oc s -1 / 2 The center-of-mass energy span corresponds to 
the RHIC energy and each curve shows the results at four different cms energies as indicated on 
the figure. 
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Fig.6 
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